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Abstract 

We perform a new calculation of the polarized next-to-leading order splitting functions, 
using the method developed by Curci, Furmanski and Petronzio. We confirm the results 
of the recent calculation by Mertig and van Neerven. 



1 Introduction 



The past few years have seen much progress in our knowledge about the nucleons' spin 
structure due to the experimental study of the spin asymmetries A^(x, Q 2 ) ~ gi(x, Q 2 )/ 
(x,Q 2 ) (N = p,n,d) in deep-inelastic scattering (DIS) with longitudinally polarized 
lepton beams and nucleon targets. Previous data on A\ by the SLAC-Yale collaboration 
[ l| have been succeeded by more accurate data from [2-4], which also cover a wider range 
in (x, Q 2 ), and results on A\ and Af have been published in || and || |7j, respectively. 

On the theoretical side, it has become possible to perform a complete and consis- 
tent study of longitudinally polarized DIS in next-to-leading order (NLO) of QCD, since 
recently results for the spin-dependent two-loop anomalous dimensions, needed for the 
NLO evolution of polarized parton distributions, have been calculated for the first time 
|| within the Operator Product Expansion (OPE). A first phenomenological NLO study 
has been presented in |§ , later followed by the analyses flO| . 

The calculation of the NLO anomalous dimensions or splitting functions is in general 
very complicated. This is true in particular for the polarized case, where the Dirac matrix 
75 and the antisymmetric tensor e^ upa enter the calculation as projectors onto definite 
helicity states of the involved particles. These (genuinely /ewr-dimensional) quantities 
lead to certain complications when dimensional regularization - which probably represents 
the only viable method of regularization in such a calculation - is used. In fact, @j was 
recently revised since an error related to the treatment of 75 was found. Although the 
results of || now fulfil a relation motivated from Supersymmetry - which appears to be 
an important constraint - it seems necessary to perform an independent calculation of the 
polarized two-loop splitting functions to check the results of ||. This is the purpose of 
this paper. 

In the unpolarized case, two different methods have been used to obtain the next- 
to-leading order splitting functions. The first calculation was performed within the 



OPE. Afterwards Curci, Furmanski and Petronzio developed a technique |i~2| , [0| which 
is as close as possible to parton model intuition since it is based explicitly on the factor- 



1 



ization properties of mass singularities in the light-like axial gauge and on the generalized 



ladder expansion fL4 |. Note that the results of |13[ fulfilled the above mentioned super- 



symmetric relation ||15|| , but were in disagreement with the first calculation |Tl| and led 
to the detection of an error in [O]. In this paper we exploit the method and results of 



T2|, pT3| to rederive the spin-dependent next-to-leading order splitting functions. To deal 



with 75 and the e tensor we use the HVBM scheme [0, which still seems to be the most 



consistent prescription |17|| . Section 2 sets the framework for the calculation, details of 
which are then given in section 3. In Section 4 we present our results. 



2 Framework 

In this section we outline the framework for our calculation. We mainly focus on the 
new features in the polarized case; more details on the method itself can be found in the 
original work (T2]. 



The general strategy consists of a rearrangement of the perturbative expansion which 
makes explicit the factorization into a part which does not contain any mass singularity 
and another one which contains all (and only) mass singularities. Fig.l represents the 
matrix element squared for polarized virtual (space-like) photon-quark scattering. The 
blob AM is expanded into 2PI kernels ACq and AKq. In the axial gauge these 2PI kernels 
have been proven Jl8[ to be finite as long as the external legs are kept unintegrated, such 
that all collinear singularities originate from the integrations over the momenta flowing in 
the lines connecting the various kernels. The generalized ladder in Fig.l can be written 
as 
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AM = AC (l + AK + AKq + ...) = AC, \ = AC AT . (1) 

ACq and Ar are then decoupled by projectors which for the longitudinally polarized 
case read (A A, AB being two polarized kernels): 

{AA)P F {AB) = (AA^iMsh) [P-P-] {^^AB k \ (2) 
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for polarized quarks and 
(AA)P G (AB) 



• ft// 1 



(3) 



for polarized gluons, putting k 2 = in the part containing the kernel AA and taking the 
pole-part (P.P.) of the projection on kernel AB. In eqs. (2), (3) i,j,k,l are spinor and 
the greek letters are Lorentz indices, n is the vector to be introduced in the axial gauge 
with n 2 = for the light-like gauge. Performing the factorization of mass singularities 



12| , which in dimensional regularization (d = 4 — 2e) appear as poles in 1/e, one ends up 



with the contribution to the (partonic) deep-inelastic structure function g^. 



Q 2 



Q 2 



AT(x,a, 



(4) 



gt(— , x, a 8 , -) = AC(— , x, a s ) 
fi^ e fi z 

where the convolution g) is defined as usual by 

(/®y)(a:)= / -f(z)g(-) . (5) 

In eq. (4) we have introduced the virtuality of the photon Q 2 , the unit of mass [i in 
dimensional regularization, the Bj0rken variable x = Q 2 /2pq and the strong coupling a s . 
Eq. (4) has a clear partonic interpretation: Ar(x, a s , -) describes the density of partons 
in the parent quark, is independent of the hard process considered and contains all the 
collinear singularities (poles in e), whereas AC(Q 2 / fi 2 ,x,a s ) is the (process-dependent) 
short-distance cross section. As was shown in [0, Ar does not depend on Q 2 , which 
is a consequence of the finiteness of the kernel AKq in the axial gauge |TB[, and allows 
for the derivation of a 'renormalization group' equation for AC with Ar related to the 
'anomalous dimension'. Thus Ar, to be convoluted with bare ('unrenormalized') parton 
densities which must cancel its 1/e poles, is equivalent to the respective Altarelli-Parisi 
kernels, e.g., 



Ar(x,a s , -) 



1 

1 - - 

e 



(6) 



for the non-singlet (NS) case. The final NLO expression for the (physical) spin-dependent 
structure function g^ then reads: 



1 Nf { 

9i(x,Q 2 ) = -Y,e 2 q lAq(x,Q 2 ) + Aq(x,Q 2 ) 

1 



+ 



a s (Q 2 
2tt 



AC q ® (Aq + Aq) + —AC g ® Ag 

iv f 



(x,Q 2 



(7) 



where Nf is the number of active flavors, and where in the full singlet case two short- 
distance cross sections AC q and AC g for scattering off incoming polarized quarks or 
gluons, respectively, exist. Here, the polarized parton distributions Ap = p^ —p^ (p = q, g) 



are to be evolved according to the polarized Altarelli-Parisi [2D| evolution equations which 

to NLO read (see, e.g., ^T|) 

d ot 

{Aq + Aq - Aq' - Af) = -± ( AP qq + AP qq ) ® ( Aq + Aq - Aq' - Aq 1 ) ,(8) 



dlnQ 2V * * 2tt 

d . s a 



d\nQ^ Aq ~ Aq) = ^( A ^" A ^)®(Ag-Ag) (9) 
for the NS quark densities and 

d ( AS \ o, / AP qq + AP qq - + AP qq , PS AP qg \ ( AS \ 
dlnQ 2 { Ag ) 2n { AP gq AP gg ) ® [ Ag ) 

in the singlet sector, where AS = J2 q {Aq+Aq) and the argument (x, Q 2 ) has been omitted 
from all parton densities. To NLO, all splitting functions in (8-10) have the perturbative 
expansion 



(10) 



( V 

APy = AJ* + -iAP£ . (11) 



The entries AP qq and AP qq ^ps start to be non-zero only beyond the leading order. For 
future reference it is convenient to introduce the NLO combinations 



APl q % s = AP l qq ±AP l qq (12) 

which according to (8,9) govern the NLO part of the evolution in the NS sector. AP qq ^ps 
is called the 'pure singlet' splitting function since it only appears in the singlet case. 



3 The Calculation 

Some representatives of the graphs to be evaluated in the calculation of the APA are 
shown in Fig. 2. We do not need to calculate the contributions from genuine two-loop 
graphs to the diagonal splitting functions AP 1 and AP gg , which are ~ 5(1 — x), since 
these are the same as for the unpolarized case. 

Before giving details of the calculation, we note that the use of the light-like (n 2 = 0) 
axial gauge in practical calculations has been a matter of debate for a long time now |22|] . 
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The great computational advantage and success it has brought for, e.g., perturbative NLO 



QCD calculations in DIS [T^, [13], ^] or jet calculus [23J has not always been matched by the 
theoretical understanding of why it worked so well [[12], 25|. The problems connected with 
the light-like axial gauge are due in the first place to the presence of spurious singularities 
in loop integrals coming from the vector propagator in this gauge, which are neither of 
ultraviolet nor of infrared origin. In (l^, [13| 'phenomenological rules' were applied to the 



problem which consisted of the subtraction of spurious poles by hand and of using the 
Cauchy principle value (CPV) prescription to deal with l/l ■ n terms (where / is some 
momentum), entailing renormalization 'constants' depending on the infinite-momentum- 
frame (IMF) variable x. In view of the complications due to the spurious poles, but 
also in view of the fact that the calculations and the prescriptions in [|l^, [l3j led to the 
correct answer, we find it very satisfactory that it turns out to be possible in our case 
to infer the effective contributions of the non-trivial p6| virtual, in particular the vertex 



correction, graphs to the polarized splitting functions from the known results |L2], |13|] for 
the unpolarized P^, such that these contributions need not be calculated all over again. 
The strategy to do this is to calculate all real emission graphs to a given process for the 
unpolarized case and to subtract their sum from the corresponding final results listed in 
|i~2| , p~3p . This difference and knowledge of the pole parts (renormalization constants) of 
the virtual corrections [O then straightforwardly allow for a determination of the result 
for each virtual correction in terms of a parametrization of the most general Dirac/Lorentz 



structure it can have |[27||, and make a transfer to the polarized case easily possible. In 



this way we avoid to have to deal with the spurious poles again but can simply rely on 



the success [28| of the approach used in the unpolarized calculation |T2|, [13 



The calculation of the real emission graphs is rather involved. This is true in particular 
for the polarized case when using the HVBM scheme since in this method the d = 4 — 2e 
dimensional space-time is explicitly decomposed into the usual four dimensions in which 
75 anticommutes with the other Dirac matrices and the — 2e dimensional part, where 
it commutes. Thus the squared matrix elements of the graphs will depend on the usual 
'd- dimensional' scalar products like l± -1% etc. (see Fig. 2a for notation), but also on \d — 4)- 
dimensional' ones, denoted by l\ ■ I2, k 2 etc. ||29||. It is most convenient to work in the 
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IMF parametrization of the momenta 



which in our case takes the form: 



V = 


( P 1 X y , P , 


n = 


( ^ - 

2P ' y ' 


k = 


/ n A; 2 + P 

xP+ - - 
\ AxP 


h = 


( ]0 Txy jz 



d -4 ) , 

r k 2 + k 2 ~ 7 \ 

% ) , (13) 

where x = kn/pn is interpreted as the IMF momentum fraction of the incoming momen- 
tum p carried by k, and k 2 = k 2 + k 2 + k 2 = k\ + k 2 is the total transverse momentum 
squared of k relative to the axis defined by p, n. We split the (d — 4)-dimensional compo- 
nents of l\ into a part parallel to those of k and a transverse part l^. According to our 
definitions, only k, l\, and I2 = p — k — l\ possess such components. When performing the 
phase space integrations one has to carefully take into account the (d — 4)-dimensional 
terms. The contribution of each graph to Ar(i,a s , =■) is given by the integration of the 
projected (see eqs.(2,3)) squared matrix elements over the phase space 

R - 1 w? xS {* - 3 / i w^ {2 ^ {1){ ' - k ~ h ~ (14) 

which is conveniently written as 
xHl - xy~ 2e [° dk 2 (-k 2 ) 1 ~ 2t f 1 dk («(1 - ~K)Y e C dw (w(l - w)y e C dv (v(l - v))^ 

J-Q 2 JO JO JO 

x( M />--)((-i- e )j( I ^M-^)(i J [ I -A.(A.(l-A.,)^)(15) 

where we have omitted trivial prefactors and defined 

k = —k 2 (l — x)kk , 

k 2 = -k 2 (l-x)k , 

ll + U = 2P(l-x)w = 2P— , 

pn 

(/?) 2 - (ID 2 = c 2 1 +v(c 2 2 -c 2 l ) = ^ 1 +ll)(l l p) , 
[II = Ai + A" (A 2 — Xi) , 
{lif = v (l - v ) ( Cl + C2 ) 2 \± (16) 
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with 



Cl,2 
Al 9 



-k 2 {\ 



X)W 



X 



(1- w){l -k)tV. 



XWK 



((/?) 2 - (llf - c lC2 ) =F (d + c^(l-«)(l-A>(l-,) . 



Note that the last two integrals in (15) are all unity if no dependence on (d— 4)-dimensional 
scalar products occurs, which of course is always the case in the unpolarized situation. 
If present, such (d — 4) -dimensional terms only give contributions proportional to e after 
the last three integrals in (15) have been performed. 



Following [0, [L3| we will regularize infrared divergencies appearing at l\n — > or 
l%n —>■ (w —>■ or w —>■ 1), which are typical of the axial gauge, by the CPV prescription 
(see above): 

1 In 

(17) 



1 

In 



(In) 2 + 5 2 (pn) 2 

All the resulting divergencies of this type can then be transformed into the basic integrals 



1 yWy 

dy 2 . x2 

o y z + cr 



0,1) 



Jo, I\ have to cancel out in the final answer which, as well as the finiteness of each 
graph in the axial gauge \T2, |TS| before the final k 2 integration is performed, helps when 
extracting the contribution of the virtual corrections by the strategy outlined above. As 
mentioned above, the renormalization constants depend on the IMF fractions x or 1 — x 
when the CPV prescription is used |12j| . We finally note that whenever considering a 
genuine ladder graph with two parallel rungs, subtraction of the 'doubly collinear' graph 



(see Fig. 2) is required within the method of \12\. The result for this is given by convoluting 
the n-dimensional leading order splitting function standing for the upper part with the 
four- dimensional one representing the lower part of the diagram, and including a factor 
(1 — x)~ e from phase space in the convolution. In 4 — 2e dimensions the polarized LO 



splitting functions read for x ^ 1 in the HVBM scheme |30 

/ 1 + x 2 x 



x) 



\1-X J 

2T R N f (2x - 1 - 2e(l - ar)) 
C F (2 - x + 2e(l - x)) , 



AP 9 ° 9 (i,e) = 2C A - 2x + 1 + 2e(l - x)) , (19) 

where CV = 4/3, = 3, Tr = 1/2 and Afy is the number of active flavors. 



4 Results 



In the normalization of [|TJ|, [13j our MS results read: 



A^sW = P^ 5 (x)-2A)CV(l-x) , (20) 

APi, PS (x) = AP^ s (x) , 

APi(x) = AP(&(x)+4C F (l-x)®AP° 9 (x) , 

AP^x) = AP#(x)-4C F (l-a;)®AP 9 Va:) , 

AP g » = AP»(x) , (21) 

where (3 = UC A /3-AT R N f /3 and APP-(x) = APP.(x,0) (see eq. (19)). The AP/i ] JHfl 
are the results of ||, and Pqq%s can t> e found in fll2"f . As was already discussed in P, [32]j 
and indicated in eq. (20), the '+' and '— ' combinations of the NS splitting functions 
as defined in (12) interchange their role in the polarized case, such that, according to 
eqs.(8,12,20), the combination AP^ NS = Pq q —Pqq — 2^ Cp(l — x) would govern the 
(Revolution of, e.g., the polarized NS quark combination 



AA 3 (x, Q 2 ) = (Aw + Aw - Ad - Ad) (x, Q 2 



Since the first moment (integral) of the latter corresponds to the nucleon matrix element 
of the NS axial vector current q^^^^X^q which is conserved, it has to be Q 2 -independent 
f33| . Keeping in mind that the first moment of the unpolarized P q l q — P*= vanishes already 
due to fermion number conservation Hl2"| , it becomes obvious that the additional term 
— 2/3 CV(l — x) in (20) spoils the Q 2 independence of the first moment of AA 3 (x, Q 2 ). 
It is therefore necessary to perform a factorization scheme transformation to the results 



in (20,21) in order to remove this additional term which, as pointed out in |30|, |32|j , 
is typical of the HVBM scheme with its not fully anticommuting 75 and trivially would 
not be present in a scheme with a fully anticommuting 75 since then the two 75 matrices 
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appearing in the relevant graphs could be removed by anticommuting them towards each 
other and using 75 = 1 (cf. Fig2a). We note that all these observations were already made 
in the original calculation of || in the OPE where, however, the removal of the additional 
term ~ (1 — x) corresponds to a finite renormalization rather than a factorization scheme 



transformation. It is nice to recover this analogy between our results and those of P, |35 



The factorization scheme transformation for AP q ^ NS also affects the singlet sector since, 
according to eq. (10), AP^+ NS +AP^ PS occurs in the evolution of the NLO quark singlet 
AS. The transformation reads in general (see, e.g., |3E[1): 

^Pqq,NS = ^Pqq,NS ~ ^Po z qq , 
^Plq,PS = ^Pqq,PS J 

AP l qg = AP l qg +Az qq ®AP G qg , 
AP l gq = AP l gq -Az qq ®AP° gq , 

APi = AP g \ , (22) 

where the APA now are the NLO splitting functions on the left-hand-sides of eqs. (20,21) 
and the APA are the new splitting functions after the scheme transformation. One im- 
mediately sees that the choice 

Zqq = ~CV(1 - X) (23) 

leads to AP q l ^ NS = P qq ^ N g and thus now yields the required vanishing of the first moment 
of AP qq + NS . Even more, the transformation (22,23) removes all additional terms on the 
right-hand-side of (20,21) simultaneously, bringing our final result into complete agree- 
ment with the revised one of M. We finally note that the above factorization scheme 
transformation also changes the quark short distance cross section (coefficient function) 
AC q in (7), since the combination 

9 A P l,± 
AC ^_2AP^ 



must be independent of the choice of the factorization scheme convention | 36fl . As was 
shown in [33, |3j], |32|, only after the transformation (22,23) takes AC q in the MS scheme 
the form 



AC q (x) = C t 



2 , \n(l-x)\ 3 1 l + x\ 

1+x 2 ) — r lnx + 

1 J \ l-x , 2(1 -x)+ 1-x 



/9 7T 2 \ 1 

+ 2 + x-(- + -U(l-x) , (24) 

i.e., becomes the previously calculated 0(a s ) quark- correct ion to g\ giving rise to, 
e.g., the correct first order correction 1 — a s /it to the Bj0rken sum rule. In eq. (24), 

[ 1 dzf(z)(g(z)) + = [ 1 dz(f(z)-f(l))g(z) . 
Jo Jo 

For completeness we note that the NLO MS gluonic short distance cross section AC g in 
eq. (7) remains unaffected by the transformation (22,23) and reads (see, e.g., 0) 



AC g (x) = 2T R N f 



(2x - 1) ( In - — - - 1 ) + 2(1 - x) 



x 



(25) 



Our final results for the polarized MS NLO splitting functions are given by: 



^^qq^S 



APi(x) 



ap;(x) 



r qq,NS 1 



C F T R N, 



AP^x) = C F T R N f 



2(1 -x) -2(l-3x) lnx-2(l + x) ln 2 x 



- 22 + 27a; - 9 In a; + 8 (1 - x) ln(l - x) 



+ -5p qg {x) (4 ln 2 (l - x) - 8 ln(l - x) lnx + 2 In 2 x - 8((2) 



-C A T R N f 



2 (12 - llx) - 8 (1 - x) ln(l - x) + 2 (1 + 8x) lnx 



-2 (ln 2 (l -x)- C(2)) Sp qg (x) - (24 - 3 In 2 x) 5p, 
C F T R N f 



qg v 



4 4 

-~(z + 4) - -Sp gg (x) ln(l -x) 



+C 2 



2 - I ( 4 - x) lnx - 5p gq (-x) ln(l - x) 



+ ( -4 - ln 2 (l - x) + - In 2 x ) <fy g(? (x) 



(4 - 13x) In x + - (10 + x) ln(l - x) + - (41 + 35x) 
3 9 



(26) 
(27) 



(2? 



(29) 



+ - (-2I X + 3\n 2 x) 5p gq (-x) + (ln 2 (l-x) -21n(l-x) lnx - ((2)) 5p gq (x) 



-C A T R N f 
-C F T R N f 



4 20 4 

4 (1 - x) + - (1 + x) mx + —5p gg (x) + -5(1 - x) 



10 (1 - x) + 2 (5 - x) lnx + 2 (1 + x) In 2 x + 5(1 - x) 
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+c 2 A 



1 19 

- (29 - 67x) Inx - — (1 - x) + 4 (1 + x) In 2 x - 2IJp gg (-x) 



67 



+ — - 4 ln(l - x) \nx + In 2 x - 2£(2) <^(x) + 3C(3) + - U(l - x) 



where, as mentioned above, the unpolarized NS pieces Ph%s can be found in Jl2| and 



Sp qg (x) 

5p gg (x) 



2x - 1 
2 - x 

1 



2x + 1 . 



(31) 



(1-x). 

Furthermore we have in eqs. (26-30) C(2) = vr 2 /6, C(3) « 1.202057 and 

Jx/(l+x) Z \ Z ) 

For relating our results to those of the relation 

4 = -2Li 2 (-z) - 2 lnxln(l + x) + - In 2 x - C(2) 

2 

is needed, where Li 2 (x) is the Dilogarithm |j4Q| . In (30), the contributions ~ 5(1 — x) to 
AP 1 are the same as those for the unpolarized P g [^TJ; they lead to satisfaction of the 



constraint 42 



1 dxAP^x) = ^ 



^C 2 A -C F T R N f - h -C A T R N } , 



valid in the MS scheme. 



In conclusion, our calculation, which was based on the approach of [0 and on using 
the HVBM |TB[ prescription for 75, has confirmed the recent results of || for the spin- 
dependent two-loop splitting functions AP^-(x). Our results also once more demonstrate 



the usefulness and applicability of the method of and the light-like axial gauge in 
perturbative QCD calculations. 
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Figure Captions 

Fig.l The matrix element squared for polarized photon-quark interaction, its expansion 
in terms of 2PI kernels ACq and AK , and its final factorized form. 

Fig2. Some representative Feynman graphs to be evaluated in the calculation of a: 
(A)Pi, b: (A)P^ (A)P^ PS , c: (A^, d: (A)P g \, e,f: (A)P g 1 g . Subtraction 
of 'doubly collinear' graphs is indicated. 
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